
function J=J_bvg(U1,U0,theta,sgi_nodes,DS,DB)

N=length(DS)/2;

J=theta(2)*U1.*(DS(N+(1:N),1)*sgi_nodes'+DB(N+(1:N),1))-...
    theta(1)*U0.*(DS(1:N,1)*sgi_nodes'+DB(1:N,1));
